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Clump mass function at an early stage of molecular cloud 
evolution: I. A statistical approach 
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ABSTRACT 

We derive the mass function of condensations (clumps) which were formed through a 
turbulent cascade over a range of spatial scales L ^ 20 pc during early, predominantly 
turbulent evolution of a molecular cloud. The approach rests upon the assumption of 
a statistical clump mass-density relationship n ex 771^ with a scale dependence of the 
exponent x obtained from equipartition relations between various forms of energy of 
clumps. The derived clump mass function (CIMF) could be represented by series of 2 
or 3 power laws, depending on the chosen equipartition relation, the velocity scaling 
index and the type of turbulent forcing. The high-mass CIMF exhibits an average slope 
r ~ — 1, typical for fractal clouds, whereas its intermediate-mass part is shallower or 
flattened, in agreement with some observational studies. 
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1 INTRODUCTION 

Understanding of the formation and the evolution of molec- 
ular clouds (MC) is of key significance in the theory of 
star formation. The early MC evolution, prior to subse- 
quent processes of active star formation, allows for simpli- 
fied physical modelling due to lack of feedback from the 
emerging s tars. Recent numerical simulat i ons shed light on 
this epoch (IVazguez-Semadeni et al.ll2007l : iHennebelle et al.l 



vation s of nearby MCs and MC co mplexes. iTesti fc Sargenj 
(| 19981 ') and I Johnstone et all (|200ll 'l derive CMFs that could 
be fitted by a single power-law function of slopes from 
-1.1 to -1.6 in the mass range m > O.SAf©, in agreement 
wit h the Salpet er high-mass slope of the IMF F — —1.3 
I Salpeted 1 19551). Other authors argue for a lognormal 
( Enoch et al.l 120081 ) or a two power-l aw shape of the CMF 
that mimics the IMF even better |Motte. Andre fc Neril 



Mi hMMleeAnlMm. Its chara cl^^I^t ic stages could be 
summarized as follows: i) convergent flows in the warm neu- 
tral medium lead to local compressions and non-linear insta- 
bilities; ii) turbulent domains (clouds) of cold molecular gas 
form in the dense regions; iii) self-gravity in the cloud takes 
slowly over and local sites of gravitational collapse emerge; 
iv) global contraction of the cloud starts (for a rev iew, see 
iMac Low fc Klessen|[2003 : IVazauez-Semadenill201(]| ').ll 

Stars begin to form at stage iii) of the MC 
evolution, as compressed cloud regions (usually la- 
belled 'clumps') evolve and fragment to (prestellar) 
cores of typical siz e 0.01 < I < 0.1 pc and densities 
n > 10-* - 10^ cm-3 iBergin fc Tafallall2007l ). Therefore the 
core mass function (CMF) is considered as a clue to the long- 
standing problem of the stellar initial mass function (IMF) 
and its variations in different environments and star forma- 
tion regions. Indeed, the correspondence between the CMF 
and the IMF is well established from dust continuum obser- 
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Nutter fc Ward-Thompson 2007 ). Their results are con- 
firmed bv lAlves. Lombardi fc Ladal (|2007l ) who used a more 
reliable approach to obtain core masses from dust extinction 
measurements toward stellar background. The CMF in the 
Pipe nebula derived by them is very similar to the IMF 
but shifted to larger masses by a factor of 4 - it appears 
that there is one-to-one correspondence between stars and 
prestellar cores, assuming star-formation efficiency (SFE) 
of ~ 25%. Thus the above-mentioned works suggest that 
the distinct dense cores are direct progenitors of stars and 
that the IMF and its characteristic mass are determined 
by turbulent fragment ation processes in MCs a n d their 
funda mental physics (jPadoan fc Nordlundl |2002| : iLarsonl 
l20 05f). The potential cave ats of this approach are discussed 
bv lciark. Klessen fc Bonn cU (200;^). 



However, the origin of the CMF from the mass distribu- 
tion of the initially formed MC clumps is an issue that still 
needs further elucidation. The clump mass function (CIMF) 
derived from CO maps exhibits significantly shallower slope 
—0.6 > r > —0.85 than the IMF and lacks a characteris- 
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tic m ass (|Blit j|l993l : iHeithausen et al.|[l998l : iKramer et al.1 
119981 ) in contrast to the characteristic mass Mch ~ 0.5Mq 
of the IMF. Later observational studies allowed for more 
detailed mapping of MCs. Emission from CO molecules 
was found to trace lower density cloud regions. The use of 
other tracers like C^^O and ^'^CO revealed structures with 
n ~ 10* cm~"^ b ut essentially larger (I ^ 0.1 — 0.5 pc ) than 
prest ellar cores (jOnishi et al.l 119961 : iTachihara et al.l l200d . 
[2003). Such compact clumps ('dense MC cores') encompass 
about 10% of the cloud mass. The derived ClMFs could be 
represented by a single power-law or a combination of 2 or 
3 power-law functions, with a steeper high-mass tail in com- 
parison to the IMF. 

Numerical simulations of clump/core formation caused 
by supersonic turbulent fragmentation yield a CIMF/CMF 
characterized by a continuum of sl opes, with a high-mass 
tail steeper than F = —1.3 (jBallcstc ros-Paredes et al.ll2006l : 
ISchmidt et al.ll2010l ). The velocity dispersion affects signif- 
icantly the shape of the time-averaged mass distribution 
- its characteristic mass decreases and the total number 
of clumps/cores grows with increasing sonic Mach num- 
ber M ( Ballcstc ros-Paredes et all [JOOa ') . The type of the 
applied turbulent forcing plays also important role. Com- 
pressive forcing produces a shallower high-mass part of the 
clump/ core mass distributi on compared to the solenoidal 
regime (jSchmidt et al.ll20ld ) . Apparently, a purely turbulent 
origin of the CIMF is inconsistent with a single power-law 
behaviour. 

This work presents a novel approach for the deriva- 
tion of the CIMF as a superposition of mass distributions 
of clumps, generated by turbulent shocks within a range 
of spatia l scales 0.5 < L < 20 pc . It is based on the 
study of iDonkov. Veltchev fc KlessenI (|201ll . hereafter Pa- 
per I) wherein clumps are defined as condensations formed 
through a turbulent cascade during the early MC evolution. 
Our starting point is the power-law relationship n oc m^ be- 
tween clump masses and densities which was substantiated 
in Paper I considering equipartitions relations between var- 
ious forms of energy. In Section [5] we list our basic physical 
assumptions, introduce the parameters of the model, and 
demonstrate how they fit the observed MC structure. The 
method to derive the CIMF is described and illustrated in 
Section [S] The results for different choices of clump energy 
equipartition are presented in Section|4l Section[5]contains a 
discussion on our model predictions in view of observational 
fits of the CIMF and some numerical results. Our conclu- 
sions are summarized in Section [6] 



iKlessen fc Hennebellell20ld ). The lower limit is close to the 
transonic scale and is set to provide generated clumps with 
sizes within the inertial range (see Sect. 4.1 in Paper I). 
Adopting a largest scale Lup ^ 20 pc, we ensure that the 
gas is mainly molecular and isothermal with typical temper- 
ature T = 10 K. 

Turbulent velocity dispersion u and mean mass density 
(p) are assu med to scale a ccording to "Larson's first and 
second laws" (|Larsonlll98ll ): 



u = 1.1 L [kms 
(p) = 13.6 X 10"^^ 



[gem 



(1) 
(2) 



where the coefficient in the second equation is derived adopt- 
ing a mean molecular mass p. = 2.4. The problem whether 
the power-law indic es a and /3 are interdep endent is still 
not resolved (see e.g. lBallesteros-Paredesll2006l l . We choose a 
fixed value of beta in the range [0.33, 0.65] whereas a = a{L) 
is derived self-similarly from the assumption of mass-density 
relationship for clumps generated at a given scale L (see 
Sect. lOI) . 

The scaling of the mean magnetic field B is obtained 
from its relation to the mean mass density. It is widely 
adopted that B oc (p)"'^ which i s verified from an exten- 
sive survey of observational data (jCrutcherll 19991 ) . Thus we 
get: 



B = 50 L" 



[pG] 



(3) 



where the value B{L — 1 pc) — 5 pG is chosen from the 
magnetic field scaling in ICrutcheiJ (|l999l . Fig. 1 there), as- 
suming a mean density scaling according to "Larson's second 
law" (equation [21) , with a = — 1. 



2.2 Clump density distribution 

Numerical simulations of supersonic turbul e nt flows (e.g. 
KlessenI [2OO0I: iLi. Klessen fc Mac Low! l2003l : iPadoan et al] 
20071 : iFederrath et al.ll20ld ) show that the volumetric distri- 
bution of (mass) density p is described statistically through 
a standard lognormal probability density function (pdf): 



p{s) da - 



exp 



5 peak 



ds 



(4) 



where s = ln(p/(p)), cr is the standard deviation (stddev) 
and Speak is the peak position: 

(5) 



cr^ = ln(l + 6^X^ 



2 STATISTICAL DESCRIPTION OF MC 
STRUCTURE 

A detailed description of the physical framework for this 
study is given in Paper I. Below we summarize it and illus- 
trate the signiflcance of the free parameters of the model. 



2.1 Scaling laws 

To describe the early evolution of MCs we consider fully 
developed turbulence creating density structures at any 
scale in the range Lup > L > 0.5 pc through a cas- 
cade possibly driven by the very process of cloud formation 



Speak 



(6) 



The dependence of o on the Mach number deter- 
mines its scaling through the choice of /3 (equation [T]). 
The turbulence forcing parameter b varies in the range 
0.2—1.0, depending on the driving type (IKritsuk et al.ll2007l : 
IFederrath. Klessen fc Schmidt|[20o3 : IFederrath et allbOld ) . 



In this paper we adopt the range of values 0.33 ^ b ^ 
0.55 which correspond to tr ansition from purely solenoidal 
to pu rely compressive case (jFederrath. Klessen fc Schmidtl 
l2008h . 

In our approach, the clumps generated at a given scale 
L axe represented by a group of statistical objects ('average 
clump ensemble') resulting from ensemble averaging over the 
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variety of Galactic clouds and cloud complexes (Paper I). 
The mass density distribution of the ensemble is lognormal, 
with a peak (equation [6l) that corresponds to its most prob- 
able member, labelled 'typical clump': 



(7) 



Pc = (p) exp(speak) = (p) exp f " Y j 

Then we define the logarithmic density range in the aver- 
age clump ensemble as ±a/2 from the most probable mass 
density p^. 



(p) exp ( Speak 



< P < (p) exp ( Speak + 



(8) 



2.3 Clump mass-density-size relationship 

The basic physical assumption about clumps is the existence 
of a mass-density relationship: 



In I -^ I = a: In 



m 
mo 



(9) 



where the power-law index x is assumed to be fixed within a 
considered clump ensemble, m is clump mass and po and mo 
are arbitrary units of normalization. Adopting the natural 
presupposition about a statistical relation between clump 
masses m, densities p and sizes /, 

(m/mo) = {p/po){l/lof , 

one obtains from equation |9] a clump size-density relation- 
ship (with size normalization unit lo) as well: 



In 



3a:: 



1 



In 



(10) 



Recalling our turbulent scenario of clump formation, the 
latter relationship should be a self-similar extension of the 
scaling of density (equation [2} , i.e. we adopt: 



3a; 



(11) 



The natural density normalization unit is the mean den- 
sity at the scale in consideration whereas the size normal- 
ization unit is chosen to be proportional to the scale size 
L: 



Po 
lo 

polo 

mo 



= (P) 
= kL 

ex exp 



(12) 
(13) 

(14) 



The last equation is derived in Paper I (equation 19 there). 
The dimensionless parameter k could be interpreted as map- 
ping resolution of the scale volume. In this work, it is taken 
to be a constant (Sect. 12. 5^ : a choice of a scale-dependent 
functional form does not yield physically meaningful solu- 
tions in terms of masses and sizes of the 'typical clump'. 



2.4 Equipartition relations 

As demonstrated in Paper I (see Section 3.1 there), the scale 
dependence of the exponent x can be derived from equipar- 
tition relations between various forms of energy per unit 



volumqj: gravitational W, kinetic (turbulent) -Ekin, ther- 
mal (internal) Eth and magnetic -Bmag- These relations yield 
a clump size-mass relationship within an 'average ensem- 
ble' at each scale L. Plotting all the ensembles on a sin- 
gle size-mass diagram, one obtains a global correlation that 
can be described by a power law and might be used as a 
diagnostic tool to verify the models through comparison 
with simulations (see Fig. 5-8 in Paper I) and observations 
(Fig. IC1|I . Best agreement with a numerical study of physics 
of clumps for med in a weakly ma gnetized turbulent medium 
with gravity (jShettv et al.ll2010l ') was achieved for choices of 
an equipartition relation that includes both gravitational 
and kinetic energy. Therefore in this work, we consider the 
following equipartition relations: 

• Equipartition of the gravitational vs. kinetic energy: 

\W\ ~ /gk^kin , (15) 

where /gk is a coefficient of proportionality. Such general 
type of equipartition is expected to hold for structures 
shaped by turbulence in which gravity gradually takes over. 
This could happen in regions where turbulence decays lo- 
cally, or where large-scale fiows accumulate material which 
eventually becomes gravitationally unstable. We adopt a 
fiducial range 1 ^ /gk ^ 4 (with an upper limit twice the 
'viria l-like' value /gk ~ 2; see e.g. IVazguez-Semadeni et al.l 
[20071) as could be expected for the early stage of the clump 
evolution. 

• Equipartition of the gravitational vs. kinetic and mag- 
netic energy: 



\W\ 



2-Ekin + En 



(16) 



This case is a form of the 'virial-like' equipartition 
\W\ ~ 2_Ekin taking into account the contribution 

of magnetic energy which is significant in some dark 

cloud cores (jCrutcher et al.l |2010| ). In fact, this equipar- 
tition is found as well through numerical simulations 
(jBallesteros-Paredes fc Vazauez-Semadenilll995l ). 

• Equipartition of the gravitational vs. kinetic and ther- 
mal energy: 



|M^|~2£kin + 2£th 



(17) 



This is also a form of the 'virial-like' case when the thermal 
component of the velocity is accounted for. The assumption 
of this equipartition is relevant at scales where clump sizes 
are typical for dense cloud cores (/ < 0.3 pc) and the internal 
energy becomes comparable to the turbulent one. 

The expressions for the forms of energy of the 'typical 
clump' are given in Appendix [X] After substituting them 
into the equipartition relations (equations I15I17]| and some 
further transformations, equations for the exponent of the 
clump mass-density relationship x are derived (Appendix 
[BJ . The obtained global size- mass correlations are illustrated 
in Appendix [C] together with observational and numerical 
reference data. The model predictions for cases |W^| ~ 2_Ekin 
(low velocity scaling index /3) and |iy| ~ 2_Ekin + 2£'th and 
\W\ ~ 4£;kin (0.42 < /? < 0.50) are in a good consistency 
with the results from the extensive observational work of 



^ Hereafter, the term 'energy' is used instead of 'energy per unit 
volume'. 
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Figure 1. Structure of MCs as traced by the exponent x{L) of the clump mass-density relationship. Observational estimates 
JLombardi. Alves &: LadajPOlOh obtained for the clouds Orion B and Taurus (dots) and California and Perseus (triangles) are com- 
pared with our model predictions, assuming equipartition relation IH^I '^ 2_Eitin (top panels) or \W\ ~ 2B]jin -|- -Emag (bottom panels) 
and for different sets of the free parameters (/3, b, k). Lines (thin: b = 0.33; thick: b = 0.55) denote the choice ft = 0.065 while shaded 
areas illustrate the effect of varying this parameter (0.02 ^ k SS 0.10) when /3 and b are fixed: (left) fi = 0.33, b = 0.33; (right) /3 = 0.42, 
b = 0.33. 



iTachihara et al.l (|2002l '). We revisit this point again in Sect.[S] 
when we comment on the results for CIMF. 



recall our self-similarity assumption: the scaling law of den- 
sity within the ensemble is an extension of the scaling law 
of the mean density (cf. equation llip . i.e. (n) oc /^^^/(i-a;) 
. Then the total mass contained within a scale L will be 



2.5 Modelled and observational structure of MCs M oc (n)L^ oc L^/t^"^). Hence, 7 = 3/(1 - x) or: 

We consider an MC as a hierarchical set of spatial scales 'y — 3 

iup ^ L > 0.5 pc which is a subset of the inertial range -^ ~ ~ 

of turbulence. Assuming that the cloud structure is condi- 
tioned by balance of energies, it can be described through 
the solutions x{L) obtained for a chosen equipartition re- 
lation. As seen from the equations in Appendix |Bl those 
solutions depend on three free parameters: velocity-scaling 
index jS, turbulent forcing parameter 6 and mapping reso- 
lution K (equation [T3| . Hence the sensitivity of x{L) to the 
choice (/3, b, k) is of critical importance. 

Observational studies of MC structure may be of help 
to restrict the para meter space that yields plausib le solu- 
tions. For instance, iLombardi. Alves fc Ladal (|20ld ) derive 
from extinction maps of several Galactic clouds and cloud 
complexes a power-law relation between the effective radius 
Rb = \/S/ti of a subregion or a set of subregions with to- 
tal area S and the mass contained in it/them: M^ oc R2. 
Since these subregions are not groups of individual clumps 
but trace the general substructure of a MC, we interprete 
their effective radii as sizes of spatial scales: L = 2Rb. 
An 'average clump ensemble' is generated at each scale L 
through a turbulent cascade in the cloud (Sect. 12.21) . We 



(18) 



Observation al estimates of x(2Rs) obta i ned in that way 
from the work of iLombardi. Alves fc Ladal (|201CI ) are com- 
pared with our model predictions x{L) in Fig. [T] Examples 
of two types of clouds are selected: with 'shallow' (Orion B, 
California) and with 'steep' (Taurus, Perseus) structure in 
terms of 2;(L). In general, best fits are obtained for lower val- 
ues of the velocity scaling index /3 and for mostly solenoidal 
forcing (6 > 0.33). The model predictions are highly sensible 
to the chosen turbulent forcing - variations of this parameter 
in the range 0.33 ^ 6 ^ 0.55 cause significant shift of the 
curve x{L) downwards and away from the zone of observa- 
tional data. On the other hand, variation of /3 affects mainly 
the shape of this curve, steepening it at larger scales. The 
value of K should be of order of 10~^ to achieve a distinction 
of substructures, significantly smaller than the spatial scale 
L and significantly larger than the scale of dissipation (Pa- 
per I). Evidently, choices of k less than a few percent yield 
solutions outside the observational range, even for 'shallow' 
clouds like Orion B. Therefore we adopt hereafter a fixed 
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Figure 2. Mass ranges of the 'average elump ensemble' vs. the scale of clump generation for different equipartition relations (columns) 
and chosen turbulent forcing (top: b = 0.55; bottom; b = 0.33) and velocity scaling index (small dots: /3 = 0.33, larger dots: /3 = 0.50). 



K = 0.065, a median mapping resolution well below 0.1. It 
allows for constructing the CIMF over a limit of confidence 
of 2-10 solar masses, depending on the assumed equiparti- 
tion. We comment on that further in Sect.U] 



3 DERIVATION OF THE CLUMP MASS 
FUNCTION 

3.1 The clump mass-scale diagram 

The turbulent parameters of each scale L: velocity dispersion 
u (equation [1} , mean density (p) (equation [2]) and forcing 
parameter b (chosen to be fixed for all scales) determine 
a lognormal density distribution of clumps pl (p) (equation 
|4]). Combining Pl{p) and the calculated x{L), one obtains a 
lognormal mass distribution of clumps pLim), 



PL{_m) 



271(7,2 



exp 



0.5 



cak — -^pcak 



Sm Sjii, peak 



/x , 



(19) 



The composite CIMF is derived as a superposition of 
the clump mass distributions pLijn), generated within the 
considered range of scales. A measure of the total number 
of clumps generated at scale L is the quantity 



A''tot {L) = — exp 



2 {x- 
a X -^ — 



l)(1^2x) 
2s2 



(20) 



as introduced in Paper I (see Eq. 17 there). It serves as 



a statistical weight for the contribution of a given scale to 
the composite CIMF. Further, the mass distributions pL{m) 
at each scale are presented through discrete sets of weights 
{NL{mj) = pL{mj)Ntot{L), j = 1, ..., n} where the range of 
masses is centred at s^, peak(i) and its limits are determined 
from the requirement: 



iVtot(L) = ^iVi(m,) 



In — - — = h„ 



(21) 



The step h,n must be chosen to be a scale-independent 
constant, in order to account correctly for the contribution of 
each scale to a given mass bin of the CIMF. That is evident 
from diagrams of scale vs. clump mass which are plotted for 
different equipartition relations and choices of (/3, b) in Fig. 
O For hm = consti, the distributions pL{rn) correspond 
to horizontal lines with a fixed linear density of dots. On 
the other hand, the choice of a constant logarithmic step 
hL ~ hi —j^ for the range of scales I/down ^ Li ^ Lup 
provides a correct counting of clumps as long as a power- 
law relationship between clump mass and scale of gener- 
ation holds. Indeed, combining equations [51 El El a^nd 1141 



one gets rric 



oc moTic 



ex {p)/gexp(-a2(3 - 22:)/22;) ex 



L"^"'"" exp(— 0-2(3 — 2x)/2a;). Transforming the index of the 
exponent in this expression by use of equations [SI and llll 
one obtains finally: 



rria oc L 



vr 



i + b^ 



1.1 



^2/3 
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where Cs — 0.186 km/s is the sound velocity at T = 10 K. 
The density scaling index a{L) is a smoothly varying func- 
tion which depends weakly on the chosen equipartition re- 
lation. As shown in Fig. [2] (left), an equipartition between 
gravitational and kinetic energy f equation I ISp yields an ap- 
proximately constant slope of the mass-scale relation. On the 
other hand, inclusion of -Emag in the equipartition (equa- 
tion I16|l causes significant change of slope (Fig. (2] right) 
at L ~ 2 pc: from a steep behaviour at small scales to a 
value of ~ 0.50, quite similar to the one derived from the 
former relation. This leads to an increased contribution of 
the small scales to the low- and intermediate-mass CIMF as 
will be demonstrated in Section (5] The mass-scale diagram 
resulting from equipartition relation between gravitational, 
kinetic and thermal energy (equation I17|l is quite similar 
(not plotted) to the one for equation 1151 

3.2 Composite CIMF 

For the derivation of the composite CIMF one must re- 
quire mass conservation throughout the entire self-similar 
structure determined by the turbulent cascade process. The 
largest scale (i.e. the whole cloud) with mass AI{Lup) is to 
conta in A/"l su bstructures of mass M{L) at given scale L < 
Lup (JElmegrecn 1997). Expressing this equation through 
the mean density (p) = M{L)/L^ and using its scaling law 
(equations [21 1111) . we get 



{p){Lup)L 

3x[L u^p ) 



1 — k(Z.up) 



Lup 

and hence 






-L^Up 



/L-^ 



(22) 



Note that Ml depends essentially on L through x{L). 
It reflects the fractal structure and can serve as a weight 
of the total number of clumps at a given scale NtotiL) for 
derivation of the composite CIMF. Then the CIMF value in 
a selected mass bin m' — Am' ^171^171' + Am' is: 



FciMF im')=Y,^Lj2 ^L [r. 



(23) 



4 RESULTS 

Brief inspection of the derived ClMFs (Fig. [3| shows that 
they can be represented by one or more power-law fits. 
By a nalogy with the fitting of the observational IMF (e.g. 
iKrou pa 2001), we define a characteristic mass M^h as the de- 
limiting value between the high-mass and the intermediate- 
mass parts of the CIMF. It is calculated as follows: i) choice 
of a trial value M^j^ ; ii) derivation of the slope (least-squares 
fit) of the high-mass CIMF with lower mass limit MJ^^^' ; iii) 
the first derived slope (i = 0) is taken as a standard to com- 
pare - if its value falls within the 3(j range of the current 
slope (i J5 1), then steps i) and ii) are being repeated. Oth- 
erwise, i.e. when the current slope differs significantly, we 
take Mch = M^^ as the characteristic mass and calculate 
the high-mass CIMF slope for all mass bins m > M^h- 

The ClMFs as derived for different equipartition re- 
lations and choices of (/3, b), with a fixed largest scale 



-Emag as found 



Lup = 20 pc, are plotted in Fig. [S] Their obvious common 
feature is the universal power-law shape of the high-mass 
part FciMF <x m^ with F ~ — 1 as expected for s elf-similar ly 
structured fractal clouds (|Elmegreenll 19971 . l2007l ) . Variations 
of r around this value are found to be minor for the equipar- 
tition between gravitational and kinetic energy (column 1 
and 3 in Fig. [3| and slightly larger (ar ~ 0.1) when other 
forms of energy are included in the equipartition. 

The mass range of the high-mass CIMF depends essen- 
tially on the assumed equipartition. In the 'virial-like case', 
without {\W\ ~ 2_Bkin) or with {\W\ ~ 2_Ekiii + 2Eth.) in- 
clusion of the thermal energy, it spans from about one to 
two orders of magnitude as the characteristic mass Mch (vi- 
olet tick mark) drops from about hundred to several so- 
lar masses with increasing the velocity scaling index /? to 
~ 0.50. If the relative weight of turbulence against grav- 
ity is increased or the magnetic energy is included in the 
equipartition, the same tendency of decreasing of M^h at 
larger /3 is evident but shifted in the range from few hun- 
dred to few tens Mq. The similarity of the characteris- 
tic mass behaviour between the results for \W\ ~ 2i?kin 
and \W\ ~ 2£;kin + 'iEth (columns 1-2 in Fig. [3]) or be- 
tween those for \W\ ~ 4_Ekin and \W\ ~ 2£'kiii + -Emag 
(columns 3-4 in Fig. [3| is remarkable. In the first case, tak- 
ing into account the contribution of the thermal energy in 
the clump energy balance leads to increase of Mch with a 
factor of 3 — 4, irrespectively of the turbulent velocity scal- 
ing. The results in the second case confirm indirectly (from 
the equivalence of the considered equipartitions) the exis- 

te nce of a statistical equipartition 2-Ekin 

bv lBallesteros-Paredes fc Vazouez-Semadenil (| 19951 ) 

In contrast to the high-mass part, the intermediate- 
mass CIMF (M < Mch) shows up a variety of shapes, corre- 
sponding to broad mass ranges from a few to hundred M© 
(Fig. [3]). When derived from the equipartitions \W\ ~ 2_Ekin 
and |VF| = 2£kin + 2-Eth, it can be fitted mainly by a single 
power-law with slopes F ~ — 0.65 as found by some authors 
for the observational CM F (|Blitj Il993l : iHeithausen et all 
119981 : iKramer et al.lll99g ). Robust fitting is possible, how- 
ever, only in cases with purely solenoidal forcing (6 — 0.33) 
and for the standard value of velocity scaling /3 = 0.33 for 
incompressible turbulence (|Kolmogorovlll94J ). Otherwise it 
fails due to two main reasons: i) the range between the lower 
mass limit of confidence (marked with green ticks in Fig. [3| 
and Mch is too small; ii) substancial part of the intermediate- 
mass range is less than 3 mass bins from the turnover of 
the mass distribution. The choice of equipartition relation 
|W| ~ 4_Bkin or \W\ ~ 2_Bkin + -Emag Icads to intermediate- 
mass ClMFs that can be described by two power-law fits. In 
the first case, the steeper part is easily recognized while fit- 
ting of the shallower one is problematic since it contains the 
turnover of the mass distribution (Fig. [3] column 3) . Their 
delimiting mass can be obtained in analogical way to the cal- 
culation of Mch. The increase of the velocity scaling index /3 
tends to equalize the slope of the steeper part to that of the 
high-mass ClMf and to steepen the shallower part. The case 
\W\ ~ 2Ekin + Emag is especially interesting because of the 
flattening of the intermediate-mass CIMF (Fig. [S] column 
4). For turbulence forcing parameter 6 = 0.33, the lowest- 
mass part exhibits even positive slopes. In view of the crite- 
ria i) and ii) introduced above, we restrain from fltting the 
intermediate-mass CIMF in that case. 
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Figure 3. Clump mass function, derived by use of different equipartition relations (columns) and for a fixed velocity scaling index /3 
(rows), choosing b = 0.33 (red) and b = 0.55 (blue). The CIMF of gravitationally unstable clumps is drawn with hatched areas. The 
lower mass limit of confidence (green tick mark) and the characteristic mass Afch (violet tick mark) are shown. Black lines denote the 
predicted slope for fractal clouds F = —1 (thick) and the intermediate-mass CIMF slope(s) for b = 0.33 (thin). 



The choice of a turbulent forcing which is a mixture of 
solenoidal and compressive mode (6 = 0.55) yields CIMFs 
of similar shape that do not depend significantly on the as- 
sumed equipartition and the velocity scaling index. The ten- 
dency with increasing /3 is toward one single power-law (i.e. 
lack of a characteristic mass) with slope F ~ — 1. Again, 
the case with presence of magnetic energy term is different, 
showing up Afch ~ 20 - 30 A/q . 

The low-mass CIMF is beyond the scope of our con- 
sideration due to the lower limit Ldown = 0.5 pc of the 
range of scales (Sect. [2TT|l and the relationship between the 
typical clump mass rric and the spatial scale of its gener- 
ation (Sect. [3T|l . At scales Ldown ^ ^ ^ 2 pc (depending 
on the chosen equipartition relation), the size of the typi- 
cal clump approaches the transonic scale which invalidates 
the assumption of supersonic turbulent fragmentation. The 
adopted velocity scaling law also fails since -Bkin ~ Eth- As 
shown in Fig. |31 the confident lower mass limit (green tick 
mark) of the derived CIMFs varies from a few to about ten 
solar masses. 



In view of the universal slope of the high-mass CIMF, 
it is worth to note the effect of varying the largest scale 
of consideration Lup- Taking the latter quantity to be the 
upper limit of the turbulent inertial range, it is propor- 
tional to the natal cloud siz e as demonstrated f rom sim- 
ulations (|Kritsuk et all l2007l : fPadoan et a!] 120071 ). Clouds 
of size Lup > 12 pc produce again CIMFs with high-mass 
slopes F ~ —1, independent on the chosen equipartition 
(Fig. 2)). Gradual decrease of Lup leads to quasi- lognormal 
shapes without a power law in the high-mass regime. On the 
other hand, intermediate-mass CIMF with its typical shal- 
lower slopes is sustained by the clump generation within a 
relatively narrow range of scales 0.5 ^ L < 5 pc which upper 
limit depends on the chosen equipartition relation. 

From the density of a chosen clump (equation[5| one can 
calculate its Jeans mass and hence estimate its gravitational 
stability. The mass functions of unstable clumps are plotted 
with hatched boxes in Fig. |31 The relative fraction of such 
objects is evidently increasing with the contribution of other 
forms of energy in the equipartition against gravity - their 
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lowest masses approach the lower mass limit of confidence 
in the case \W\ ~ 4£;kin while for \W\ ~ 2_Bkin + -Emag 
all clumps are unstable. Apparently, the mass functions of 
unstable clumps are not significantly affected by variations 
of the turbulent forcing parameter b. On the other hand, the 
fraction of unstable clumps decreases - as expected, - when 
/3 is increased. 



5 DISCUSSION 

It seems not surprising to us to obtain a slope F = — 1 for 
the high-mass CI MF, like the typical one for fractal clouds 
lElmegreenll 19971 1. in case of equipartition relation between 
the gravitational and the kinetic energy. In view of the den- 
sity and velocity scaling laws (Sect. [2TT]), turbulence tends 
to dominate against gravity at larger scales (see Appendix 
IB1|) and determines the fractal structure of the cloud and, 
hence, the high-mass CIMF. The effect on F of including the 
thermal energy in the equipartition relation is small since 
the relative weight of ^th decreases with increasing L (see 
Appendix IB3|) . Inclusion of magnetic energy in the clump 
energy balance introduces larger, although not essential un- 
certainties of the high-mass CIMF slope. This could be ex- 
plained with the strong scale dependence of the magnetic 
field (equation[3| which results in a steeper relation between 
the magnetic and kinetic energy terms (see Appendix IB2|) . 
The effect can be seen in Fig. [2] (right) - the presence of mag- 
netic field leads to a slow increase of the clump masses with 
L at small scales while at large scales the clump mass-scale 
relationship is determined mainly by the kinetic energy. 

The obtained slopes of the intermediate-mass 
ClMFs deserve special attention because of their 

i dis)similariti e s to s o me observational results on the CMF 
Onishi et al. 199^: Kramer et al. 19981: Tachihara et al.l 



l2002l : iKainulainen et al.l l201ll ). iKainulainen et all (120111 ) 
derived a single power-law CMF from extinction maps of 
nearby MCs and obtained a slope F = —0.4 ± 0.2 which 
is too shallow but possibly affected by clump blending. 
iKramer etlj] (| 19981 ') studied several MCs in hues of "CO 
and C^*0 and derived single power-law CMFs of slope 
—0.6 > F > —0.8 spanning about 2 orders of magnitude, 
over a large variety of mass ranges (from small cores up to 
clouds). Single power-law fits of similar slope result from 
equipartitions \W\ ~ 2Euin and \W\ ~ 2Skin + 2£;th (Fig. [3 
left) and - in case of larger b and/or /3, - from |Vy| ~ 4i5kin 
(Fig. [S] third column). We note that for Orion B region the 
choice \W\ ~ 2Ski„ (/3 = 0.33, b = 0.33) leads not only to 
a good agreement between the CIMF and the CMF but 
als o between the predicted MC stru cture and that derived 
by iLombardi. Alves fc Ladal (|2010l ) (Fig. [T] left top). It 
should be pointed out that the CIMF mass ranges are 
severely restricted by the lower mass limit of confidence. 
In some other cases, the intermediate-mass CIMF is ap- 
parently a combination of two power-laws (Fig. [S] right). 
The steeper part of it, when derived from equipartition 
\W\ ~ 4_Ekin, for shallow velocity scaling (0.33 sj /3 sC 0.42) 
and purely solenoidal forcing ( b — 0.33) is ag a in in general 
agreement with the results of iKramer et al.l (| 19981 ) while 
the lower-mass part flattens. Similar behaviour is found 
by iTachihara et al.l (12003) from extensive statistics of 
nearby star-forming regions although these authors obtain 



essentially steeper slope (F — —1.5) for m > 10 Mq. The 
case \W\ ~ 2_Bkin + £mag yields a combination of positive 
and flat slopes which could be considered also as a single 
flat intermediate-mass CIMF (Fig. (5] bottom right; note 
the restrictions of the lower m ass limit of c o nfiden ce) . The 
latter fitting is suggested by I Onishi et al.l (|l996l ) for the 
CMF in Taurus MC (F = 0.1). It agrees qualitatively with 
the intermediate-mass CIMF slopes derived by choosing 
values of /3 and b that fit best the observed cloud structure 
for scales 2 < L < 8 pc (Fig.[T] right bottom). A comparison 
with Fig. [2] (right bottom) shows that the corresponding 
mass ranges are 20 < tti < 100 Mp) and overlap part ially 
with the mass range of the study of lOnishi et al.l (| 19961 ). To 
sum up, the best agreement between ClMFs and observa- 
tional CMFs is apparently achieved in cases \W\ ~ 2iJkin 
(0.33 s; /3 s; 0.42) and \W\ ~ 2£;kin + 2Eth and \W\ ~ 4£kin 
(0.42 ^ /3 ^ 0.50) which is confirmed also by comparison on 
size-mass diagrams (Fig. ICip. 

Since IKramer et all (|l998h derived CMFs of single 
clouds, we may speculate that a result from a statistical 
approach, combining data from more clouds, could be a 
tw o- or three pow e r-law CMF. In that aspect, the work 
of ITachihara et al.l (|2002f ) is instructive although we still 
lack an explanation of the steep slopes they derived for 
m > 10 Mq. Their study encompasses objects of sizes, 
corresponding to our clumps (0.1 < I < 0.4 pc), but 
the variety of their dynamical state is huge: starless, star- 
forming and cluster-forming cores. In a forthcoming pa- 
per, we will address the problem, comparing the predic- 
tions of our model for chosen equipartition relation with 
CMFs derived from molecular-line and dust-continuum data 
for several Galactic MCs, refiecting various physical condi- 
tions. Within this work, we simply point out that differ- 
ent models of cloud stru c ture consistent with the study of 
iLombardi. Alves fc Ladal (|2010l ) (Fig.[Tl cf. also Fig. 1 in Pa- 
per I) lead to ClMFs which are in a good general agreement 
with some observational CMFs (Fig. [5]|. 

We obtain quasi-lognormal shapes of the high-mass 
CIMF only by restricting the upper scale Lup of the hi- 
erarchy of clumps to very small values (Fig. Q. That 
result seems to contradict the mass distributions of 
cores/clumps obtained from numerical simulations of tur- 
bulen t fragmentation for various diff erent rms Mach num- 
bers jBallesteros-Paredes et al. 2006h and choices of turbu- 



lent forcing ( Schmidt et al.ll20ld ). However, self-gravity is 



not included in those simulations while it is a main factor in 
our approach to derive CIMF from energy equipartition rela- 
tions. As already commented above, a single power-law high- 
mass CIMF is a universal feature, resulting from the balance 
between the gravitational and (mainly) turbulent energy in 
clumps. Core mass distributions of logno rmal shapes and 
with large widths were derived as well bv lDib et al.l (|2008l ) 
from simulations, that include self-gravity. We believe their 
results are different because their objects are more compact 
and probably more evolved in comparison with the clumps 
in our consideration. 

Our discussion does not include cases with high turbu- 
lent forcing parameter (1 > 6 > 0.55) where the compressive 
mode provides the mai n contribution to the turbulent energy 
(JFederrath et al.ll20ld ). Our approach is limited in that as- 
pect since the increase of the stddev of the density distri- 
bution (equation [5]) often yields unrealistically high typical 
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Figure 4. Effect of varying of the upper scale limit Lup = 12, 7, 5, 2 pc (decreasing line-width) on the clump mass function, derived 
by use of different equipartition relations (columns) and choosing b = 0.33 (bottom) and b = 0.55 (top), whereas /9 = 0.42 is fixed. The 
slope r = — 1 (thick black line) and the lower mass limit of confidence (green tick mark) are plotted for comparison. 



clump masses rric in relation to the mass included within the 
volume of a given scale L^ . That constrains the spatial range 
for a plausible derivation of the CIMF. Such solutions are 
obtained only for low /3 ~ 0.33 and do not differ significantly 
from those derived for b = 0.55 (Fig. O upper panels): sin- 
gle power-low with slope about —1. Nevertheless, the studied 
cases for h = 0.55 are representative enough for the effect of 
the compressive mode on the CIMF. 



6 SUMMARY 

By use of a statistical approach, we derived mass functions 
of condensations (clumps) which were formed through a tur- 
bulent cascade over a range of spatial scales L ^ 20 pc dur- 
ing the early MC evolution. Clumps are considered within 
the framework of Paper I: as ensembles of objects in a state 
of equipartition between gravity and other forms of energy 
and obeying a power-law mass-density relationship n oc m^ . 
The functional form x = x{L) is determined by the cho- 
sen equipartition relation and the free parameters of the 
model: velocity scaling index 0.33 ^ /3 ^ 0.65 and turbulent 
forcing parameter 0.33 ^ 6 ^ 0.55. The clump mass distri- 
bution at a fixed scale was obtained from the assumed log- 
normal density distribution and then the composite clump 
mass function (CIMF) was derived by superposition of the 
clump mass distributions generated at the various different 
scales, assuming self-similar cloud structure. 

The obtained ClMFs for different equipartition relations 
could be represented by series of power-law functions as 
intermediate-mass and high-mass parts are distinguished, 



with a characteristic mass Mch that varies from a few to 
a few hundred Mq. The high-mass CIMF can be fitted by 
a powe r-law of average slope F ~ — 1, typical for fractal 
clouds (|Elmegreenl 119971 ). with some variations (or ~ 0.1) 
when magnetic energy is included in the energy balance of 
the clumps. When derived from a 'virial-like' equipartition 
without {\W\ ~ 2_Bkin) or with accounting for the ther- 
mal component of the velocity {\W\ ~ 2_Bkin + 2i?th), the 
intermediate-mass CIMF could be represented by a single 
power-law of slope F ~ —0.65, in agreement with some ob- 
servational clump mass functions (CMFs). Increase of the 
contribution of turbulent {\W\ ~ 4i5kin) or magnetic en- 
ergy {\W\ ~ 2ii'i:in + -Emag) agaiust gravity in the clump 
energy balance leads to an intermediate-mass CIMF which 
is a combination of two power-laws, except in the case of 
large velocity scaling index (/3 > 0.50). The slope of the 
steeper part varies in a narrow range —0.7 > F > —0.9 
depending on the adopted equipartition. The other power 
law tends to flatten in case of purely solenoidal turbulent 
forcing (b = 0.33) and even has a positive slope when the 
equipartition | W| ~ 2£'kin + Smag is adopted. 

Careful comparison with observationally derived mass 
functions of clumps, considerably larger than low-mass 
prestellar cores, would demonstrate the ability of our model 
to match the variety of physical conditions in Galactic star- 
forming clouds. 

Acknowledgement: T.V. acknowledges support by the 
Deutsche Forschungsgemeinschaft (DFG) under grant KL 
1358/13-1. 
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APPENDIX A: CLUMP ENERGIES 

Below we list the formulae for the terms of energies of the 
'typical clump' (subscript 'c') used in this work. 

• Gravitational energy: 

|h^N.4g£^p. (A1) 

The coefHcient z^ accounts for the contribution of the mass 
outside the clump, typically varying between 1 (vanishing 
gravitational influence) an d 2 (strong gravitational influenc e 
from the external cloud) (jBallesteros-Paredes et al.l 120091 ). 
That range is applicable for scales below the sizes of giant 
MC and we adopt Zc = 1.5 (moderate external gravitational 
influence) in all considered cases. 

• Kinetic (turbulent) energy: 



P 1 2 1 2 

iikin = -^PcUc = -pcUo 



Ic 

1 pc 



2/3 



(A2) 



where the typical clump velocity scales like the rms velocity. 
Mo = 1.1 km/s (equation [ij. 

• Internal (thermal) energy: 

3 V 
Eti, = ^-pcT , (A3) 

where 5R is the gas constant. 

• Magnetic energy: 

£^ag-^^^, (A4) 

where the magnetic field B scales according to equation [3l 



Bl Gravitational vs. kinetic energy 

From equation 1151 
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B2 Gravitational vs. kinetic and magnetic energy 

From equation 1161 
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B3 Gravitational vs. kinetic and thermal energy 

From equation 1171 

' 4a; -f 2 a: - 1 
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APPENDIX C: SIZE-MASS DIAGRAMS 

Plotted on size-mass diagrams, the 'average clump ensem- 
bles' generated at each scale L exhibit global power-law 
correlations with slopes 7giob. In Fig. ICll they are juxta- 
posed with reference d ata from observations of cloud cores 
(jTachihara et al.l l2002h and simulations of cloud clumps 
for med in a weakly mag netized turbulent medium with grav- 
ity (|Shettv et al.l [2010|). In the latter case, clumps have 
been delinated analogically to an observational study of MC 
structures, imposing a cut-off level on the column density 
map. 



APPENDIX B: EQUIPARTITION FUNCTIONS 

These functions are used to derive the clump mass-density 
exponent a; at a given scale. They are obtained from the 
equipartition relations (equations I15I17[I by substitutions 
from the expressions for different clump energies (Appendix 
[A]! wherein Ic and rric are excluded by use of equations [10] 
and [9] and the mass normalization unit mo is expressed from 
equation 1 141 
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Figure CI. Size-mass diagrams, composed for different equipartition relations (columns) and for a fixed velocity scaling index /3 (rows), 
choosing b = 0.33 (red) and b = 0.55 ( blue). The derived slop e for each chosen method and parameter set is given (solid black line). 
The region o f the cloud cores stu died bv lTachihara et al.l 1(2003) (green hatched areas) and the slope —1.95 obtained for clumps from the 
simulation of lShettv et al.l l|2010l ) (violet line) are plotted for comparison. 



